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Abstract. We provide some considerations on the excitation of black hole quasinormal modes 
(QNMs) in different physical scenarios. Considering a simple model in which a stream of particles 
accretes onto a black hole, we show that resonant QNM excitation by hyperaccretion requires a 
significant amount of fine-tuning, and is quite unlikely to occur in nature. Then we summarize and 
discuss present estimates of black hole QNM excitation from gravitational collapse, distorted black 
holes and head-on black hole collisions. We emphasize the areas that, in our opinion, are in urgent 
need of further investigation from the point of view of gravitational wave source modeling. 

A leading candidate source of detectable waves for Earth-based interferometers such 
as LIGO, Virgo, GEO600 and TAMA, as well as for the space-based interferometer 
LISA, is the inspiral and merger of binary black holes. The waveform should comprise 
three parts, usually referred to as inspiral, merger and ringdown. The inspiral waveform, 
originating from that part of the decaying orbit leading up to the innermost stable orbit, 
can be analyzed using post-Newtonian theory. Extensive studies of the detectability of 
this phase of the signal have been carried out, both for Earth-based 1 1 ] and for space- 
based interferometers 10] • The nature of the merger waveform is largely unknown at 
present, and is the subject of work in numerical relativity. 

The ringdown waveform originates from the distorted final black hole, and consists of 
a superposition of quasi-normal modes (QNMs). Each mode has an oscillation frequency 
and a damping time that are uniquely determined by the mass M and specific angular 
momentum j = J/M^ of the black hole [3]. The amplitudes and phases of the various 
modes are determined by the specific process that formed the final hole. 

The uniqueness of the modes' frequencies and damping times is directly related to the 
"no hair" theorem of general relativistic black holes. Therefore a reliable detection and 
accurate identification of QNMs could provide the "smoking gun" for black holes and 
an important test of general relativity in the strong-field regime [4]. 

The detectability of ringdown waves, and the accuracy with which we can measure 
their frequencies and damping times to test the no-hair theorem, depend mainly on the 
energy carried by each mode [5]. In turn, the energy distribution depends on the details 
of the merger process. Given our poor understanding of the merger phase, we have at 
best sketchy information concerning the energy distribution between different modes. 

In the first part of this paper we show by a simple toy model that QNM excitation 
by infalling matter requires a significant amount of fine-tuning, and is quite unlikely in 



astrophysically realistic scenarios. In the second part we briefly review present estimates 
of QNM excitation in various physical situations, including gravitational collapse, sim- 
ulations of single distorted black holes and head-on collisions of two black holes. This 
second part is an extended version of Sec. VB in [5]. Most of our considerations can be 
applied to the solar-mass black holes detectable by Earth-based interferometers, but the 
main motivation for this short review comes from our study of the massive black holes 
detectable by LISA |,5J. Throughout this paper we use geometrical units {G — c— 1). 

EXCITATION BY INFALLING MATTER 

Various authors |0, 0] suggested that quasinormal ringing could be resonantly excited 
by clumps of matter falling into the black hole at some appropriate rate. This rate should 
be such that the spatial separation of the clumps equals a multiple of the typical QNM 
wavelength. Here we show by a simple model that, in principle, the modes of a black 
hole can indeed be excited in this way. We also show that a simple addition of damped 
sinusoids provides a good fit of the resulting gravitational waveforms: in other words, 
we can interpret the resonant excitation of the modes in terms of interference between 
gravitational waves. For simplicity we consider clumps of mass /i much smaller than the 
black hole mass, -C M, so that we can apply a perturbative analysis. 

Consider first one single clump falling from infinity. This process was first analyzed in 
a classical paper by Davis, Ruffini, Press and Price 1 8] (hereafter DRPP). They found that 
the total radiated energy is ~ 0.0104/i^/M, most of it (0.0091/i^/M) being emitted as 
quadrupolar = 2) waves, and that the total energy spectrum is peaked at (O ~ 0.32/M, 
very close to the fundamental £ = 2 QNM frequency ft) ~ 03131 /M. 

We can now superpose the waveform and its Fourier transform for the one-particle 
case to represent two or more particles falling into a black hole and to study interference 
phenomena (we are assuming that the clumps are non-interacting, so that there is no extra 
contribution to the total energy-momentum tensor). A similar sum-over-point-particles 
approach has been used many times in the past [|9|], but the particular process we consider 
here has not been studied before. 

Let us first suppose that we drop two particles with a temporal separation of T (which, 
bearing in mind that the retarded time is the measured quantity, also means a spatial 
separation of T between the two bodies). If ^i(?,r) represents the waveform for the 
first body (normalized by /i, where /i is the mass of the infalling body), then the total 
waveform will be 

^(?,r)=M^i(f,r)+M^i(?-r,r). (1) 
Likewise the Fourier transform ^^(ft), r) will be 

^ = Ai*i+Ai«'''"^'i'i- (2) 

For A'^ bodies dropped regularly, such that the temporal difference between the dropping 
of the first and of the last is T , we have 

^(f,r)= £'M^i(r--^j,rV »I'(co,r)= ^V^^^^'^i- 0) 
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FIGURE 1. Characteristic interference pattern in the total energy radiated by two (left) and twenty 
(right) particles falling into a Schwarzschild black hole, as a function of the particles' initial separation. As 
r ^ the particles behave as a single particle of mass {Njj.), and the total energy scales like (Nn)^. The 
red, dashed horizontal line corresponds to the sum of the single-particle energies (E ^ N x 0.009 
for particles). In the two-particle case we also overplot the prediction of the two-mode ringdown 
waveform, Eq. (|8}. The blue, dot-dashed line corresponds to COr — 0.3737/M (the fundamental QNM 
frequency with £ = 2). To obtain the green, dashed line we choose the "phenomenological" value cOr = 
0.32/M (corresponding to the peak of the one-particle energy spectrum). 



The infall of a finite-sized body can also be studied with this formalism ^ . Incidentally, 
one can prove that \lJ-T!j=Q e'^-']^ < n^N^, so the total energy Ej-ad ~ radiated 
by a stream of plunging particles is always less than or equal to the energy radiated by a 
single point particle with total mass ^N. 

The total energy radiated in £ = 2 as a function of separation for two and twenty 
particles is shown in Fig. [T] Let us first consider the two-particle case (left panel). 
As we increase the distance between the two bodies the total energy decreases. At 
the first minimum (T ~ 9.9M) the total energy radiated is ~ 0.0085/i^/A/. The next 
maximum is attained at T ~ 19. 8M, the total radiated energy being ~ 0.0206/i^/M. 
Since the one-particle spectrum is highly peaked at CO ~ 0.32/M, it is tempting to explain 
this behavior as an interference phenomenon. The typical dominant wavelength of the 
fundamental QNM is A = 2;r/(y ~ 19. 6M: this means that, for two particles, destructive 
interference should occur at T ~ 9.8M and constructive interference at T ~ 19.6M, in 
good agreement with our numerical results. In the limit of large distances, interference 
becomes negligible and the total radiated energy tends to the sum of the individual 



We can think of a finite body as being composed of many particles, the separation between the particles 
being infinitesimal. The total waveform ^ will be a superposition of single-particle waveforms ^i, of 
the form 1'(f,r) — limAf^oo £^Jq^ (f — ■ The Fourier transform of this waveform is 4^ = 

iir 1 11 i(OT ■ ~ 

UmA^^oo L/=o ^^"^ I'l • By taking this limit we obtain the gravitational radiation generated by a massive 
thinbody of length T and total mass jj.. Summation of the series yields limN^ocYJJ-o ^^"^^ = ^-^wr~ ' 
sothafP^AlTjr^*!- 
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FIGURE 2. Two-particle waveforms (solid black lines) are compared with the single-particle wave- 
forms (dashed red lines). In the left panel T = 9.9M, and the ringdown is suppressed by destructive 
interference. In the right panel T — 19. 8M, and the ringdown is enhanced by constructive interference. 



particle energies (~ 0.0182/i^/M for our two particles). This interpretation of the result 
is confirmed by an explicit calculation of the gravitational waveforms (Fig.|2l). 

The situation for more than two particles is similar. Notice that in the two-particle 
case, for T < 6M there is always an enhancement on the total radiated energy (as 
compared to the sum of the single-particle energies). However, as we increase the 
number of particles, this enhancement interval decreases. So we should expect it will 
be hard to resonantly excite the black hole QNMs using a large number of particles. 

Interpreting the results in terms of ringdown waveforms - The DRPP signal is 
dominated by quasinormal ringing, so it makes sense to try and interpret the preceding 
results in terms of ringdown waveforms. To start with, we assume our signal has the 
form 

\if = @{t)Ae^'^''^in(Ort, (4) 

where @{t) is the Heaviside function. For the fundamental mode of a Schwarzschild 
black hole (Or ^ 03131 /M and (0/ ~ 0.0890/M. The amplitude A can be obtained 
imposing that the total radiated energy in £ = 2 is equal to 0.0091/i^/M. The total energy 
radiated in the £-th multipole is 



6An{l-2)\J-oo \dt 
which for the waveform ^ gives 



A\l + l)\co^ 



Setting 1 = 2 and equating this to 0.0091/i^/M gives A ~ 0.441/i, in qualitative agree- 
ment with the DRPP waveform. Let us now consider two particles by superposing two 



ringdown waveforms: 



V/ = &{t)Ae-'^'' sin cOrt + Q{t - r)Ae"®'('"^) sin (Or{t - T) . (7) 
The total energy radiated would be 

" mn{l -2)1(0, ^ ^ ^'^""^ " ^' '''' ■ 

A plot of this function is shown in the left panel of Fig. [U It is in good qualitative 
agreement with the full numerical calculation. If we use a "phenomenological" QNM 
frequency ft),- = 0.32/M (corresponding to the peak of the one-particle energy spectrum, 
and yielding an amplitude A ^ 0.515/i) instead of the "true" fundamental i = 2 QNM 
frequency (Or = 0.3737/A/, the agreement between Eq. (HI) and the numerical results 
becomes even better. 

This agreement suggests the possibility to extend our study to Kerr black holes, by 
simply assuming that (Or and cOj axe functions of the black hole's angular momentum. 
Using ® and ® we get 

^2ptcles y X 

-J^ — = 2e-^"' ( e"'^ + cos (OrT - — sin (D,r ) ~ 2( 1 + cos (OrT) . (9) 

£rad V 0)r J 

In the last step we considered rapidly rotating (near-extremal) black holes and modes 
with i = m, for which ft), — > and (Or m/{2M). Therefore, for corotating modes of 
rapidly rotating black holes the resonance is almost perfect if we choose T = 2TzM/m 
(the ratio is then 4, the maximum allowed) while destructive interference results if 
T = TzM/m. 

In summary: constructive resonance, and consequently QNM excitation, is possible 
(in principle) when a stream of particles accretes onto a black hole. Unfortunately, 
the separation between each particle required for an enhancement in the total energy 
becomes very small as the number of particles increases. Unless some delicate fine 
tuning is at work O], the separation between any two infalling particles or clumps 
should be driven by a random process, which means that on average there should be no 
resonance at all. This is confirmed by simulations of generic configurations of infalling 
matter, where quasinormal ringing is usually absent from the waveforms [.10.1 . 



ENERGY DISTRIBUTION IN DIFFERENT SCENARIOS 

In this Section we present a brief review of results on QNM excitation in different 
physical scenarios, including gravitational collapse, evolutions of single distorted black 
holes and head-on black hole collisions. Far from being complete or comprehensive, our 
review aims to emphasize areas that require further investigation from the point of view 
of gravitational wave source modeling. 

Gravitational collapse - There are just a few calculations of gravitational wave emis- 
sion from rotating gravitational collapse to a black hole. Perturbative calculations were 



first carried out by Cunningham, Price and Moncrief and later improved 

upon by Seidel and collaborators [ 14]. These studies suggest that gravitational waves are 
mainly generated in the region where the Zerilli (or Regge- Wheeler) potential is large, 
and that the signal is always dominated by quasinormal ringing of the finally formed 
black hole. Simplified simulations based on a free-fall (Oppenheimer-Snyder) collapse 
model yield a rather pessimistic energy output, with a typical core-collapse radiating 
up to £■ ~ 10~^M in gravitational wave energy, and most of the radiation being emitted 
in modes with i = 2. Radiation in £ = 3 is typically two to three orders of magnitude 
smaller than the i = 2 radiation (see eg. Fig. 9 in fllj). Cunningham et al. also provide 
analytical expressions for the energy radiated in different multipoles as a function of the 
initial quadrupole deformation of the star and of the initial collapse radius [ fTlL [l2|l. 

For a long time the only nonperturbative, axisymmetric calculation of gravitational 
wave emission from collapse has been provided by the seminal work of Stark and Piran 
111 51 . Surprisingly, the waveform resembles that emitted by a point particle falling into 
a black hole, but with a reduced amplitude (we will see below that a similar result 
emerges from head-on black hole collisions in full general relativity). As expected for 
axisymmetric configurations, the amplitude of the cross polarization mode is always < 
0.2 that of the even mode, and the emission in the plus mode is dominantly quadrupolar, 
with a sin 9 angular dependence (the cross polarization has a dependence of the form 
cos Q sin^ 0). The total energy emitted is quite low; it increases with rotation rate, ranging 
from ~ lO^^M for j = to ~ 7 x IQ^^M as j ^ 1 . Rotational effects halt the collapse for 
some critical value of j = which is very close to unity, and depends on the amount 
of artificial pressure reduction used to trigger the collapse. Stark and Piran find that the 
energy emitted scales as for low j: E/M ~ min{1.4 x lO^^j^, Cmax} for < j < jcm, 
where Cmax ~ 10^^ again depends on the amount of artificial pressure reduction used 
to trigger the collapse. The /i+ waveform is very well fitted by a combination of the 
two lowest QNMs with £ — 2 and it is only weakly sensitive to j, reflecting the weak 
dependence of the slowly-damped Kerr QNM frequencies on j. 

This two-dimensional calculation has recently been improved upon using a three- 
dimensional code [16], but still keeping a high degree of axisymmetry (so only modes 
with m = are excited). Ref. [ 16] picks as the initial configuration the most rapidly rotat- 
ing, dynamically unstable model described by a polytropic EOS with F = 2 and K = 100, 
having a dimensionless rotation rate ~ 0.54, and triggers collapse reducing the pressure 
by less than about 2%. Gravitational waves are extracted by a gauge-invariant approach 
in which the spacetime is matched with nonspherical perturbations of a Schwarzschild 
black hole. The "plus" polarization is essentially a superposition of modes with i = 2 
and i = 4, and the "cross" polarization is a superposition of modes with i = 3 and £ = 5. 
Fig. 1 in 1 16] shows that the amplitude of the i = 2 mode is about an order of magnitude 
larger than the amplitude of the i = 4 mode, and Fig. 3 shows that the cross polariza- 
tion is roughly suppressed by one order of magnitude with respect to the plus polariza- 
tion, with maximum amplitudes in a ratio |(r/M)/zx |max/|('"/^)^+lmax — 0.06 (here r 
denotes the distance between source and detector). Odd-parity perturbations should be 
zero in spacetimes with axial and equatorial symmetries: they are nonzero because of the 
rotationally induced coupling with even-parity perturbations The max- 

imum amplitude from these three-dimensional simulations is |(''/A/)/z^|jj^ax 

- 0.00225, 

about one order of magnitude smaller than the amplitude found by Stark and Piran. Cor- 



respondingly, the total energy lost to gravitational waves is £ ^ 1.45 x 1O^^(M/M0). 
This is about two orders of magnitude smaller than the estimate by Stark and Piran for 
the same value of the angular momentum, but larger than the energy losses found in 
recent calculations of rotating stellar core collapse to protoneutron stars iHtIi . 

Distorted black holes - Distorted black hole simulations shed some light on the 
(crucial) issue of the dependence of QNM excitation on the initial data: that is, on the 
"shape" and "size" of the black hole distortion. 

Ref. [18] carries out a pioneering study of linearized gravitational perturbations of the 
Kerr spacetime. Given initial data with a certain angular dependence - that is, assigned 
values of (/, m) - their study finds that modes with positive and negative m are always 
excited together. The simultaneous excitation of corotating and counterrotating modes 
depends on a certain reflection- symmetry property of the Kerr QNM spectrum [5], and is 
confirmed by more recent time-evolutions of scalar perturbations of the Kerr spacetime 
|[l9]. In addition, these new simulations show that the relative excitation of corotating 
and counterrotating modes depends on the radial profile of the initial data, confirming 
expectations based on analytic investigations of the same problem [2flJ. Unfortunately, 
detailed studies of the excitation of gravitational perturbations of Kerr black holes in the 
perturbative regime are still lacking. 

Ref. [21] provides one of the first attempts to extract gravitational waveforms from 
non-linear simulations of distorted black holes. They use a two-dimensional code to 
evolve initial data corresponding to a single, nonrotating black hole superimposed with 
time-symmetric gravitational waves (Brill waves), extract gravitational waves using 
gauge-invariant perturbation theory and show that, for low-amplitude Brill waves, non- 
perturbative and perturbative waveforms with £ = 2 and £ = 4 are in good agreement. 

A more complete analysis is due to Anninos et al. [22]. They use three different initial 
data sets, corresponding to very different physical situations: non-rotating black holes 
distorted by time- symmetric Brill waves, distorted rotating black holes, and the time- 
symmetric "two black hole" Misner data. In all three cases the system quickly settles 
down to a nearly spherical (or oblate, for rotating holes) configuration. They study 
the dynamics of the black hole horizon monitoring the ratio of polar and equatorial 
circumference Cr = Cp/Ce, and find that the horizon oscillates at the QNM frequencies 
of the final black hole. Oscillations with i = 2 and i = 4 are well fitted by a superposition 
of the fundamental QNM and the first overtone. In the "black hole plus (small) Brill 
wave" case, even when the initial data contain a significant component with i = 4 (that 
is, when the index n in the analytic expression of the Brill wave is equal to 4) the i = 4 
component of the horizon distortion is a factor ~ 10^ smaller than the i = 2 component 
(Figs. 2 and 3 in lE2ll ). For Misner initial data, corresponding to colliding black holes, 
£ = 2 waves are also strongly dominant. An important conclusion of Ref. ll22ll is that the 
dynamics of the apparent horizon geometry can be used not only to find the fundamental 
QNM frequency of the hole, but also its mass and angular momentum. 



Brandt and Seidel 112311 present an extensive study of distorted, rotating "Kerr-like" 
black holes with a wide range of rotation parameters (as high as j ~ 0.87). Quasi- 
normal ringing is observed both in the even and in the odd components of the emit- 
ted radiation. The angular momentum parameter j is well approximated by the simple 
fit j — — ( — 1.55 -1-2. 55Cr)^ (the fit being generally accurate within ~ 2.5%): this 



means that a measurement of the horizon geometry provides a value of the rotation pa- 
rameter. For black holes with j > \/3/2, a Kerr black hole cannot be embedded in flat 
space: an alternative method to obtain j consists in measuring the angle at which the 
embedding ceases to exist. 

A crucial result of Ref. ll23ll is that the black hole dynamics depend critically on the 
properties of the initial data sets. For example, initial data corresponding to nonlinear, 
odd perturbations produce oscillations of the Gaussian curvature at frequencies given by 
linear combinations of QNMs with £ = 3 and i = 5, because of beating phenomena. In 
Sec. rVB, they present a method to "dig out" frequencies with stronger damping from 
the slowly-damped modes in a Fourier transform, which can find useful applications in 
gravitational-wave data analysis. Sec. V of [23] presents waveforms corresponding to 
different initial distortions of the hole. Odd-parity distortions of Schwarzschild black 
holes (run "ol" in their terminology) produce a waveform in which 99.9% of the energy 
is radiated in £ = 3, 0.05% of the radiation goes into i = 2, 0.02% into £ = 4 and 0.0001 % 
into i = 5; but for different "shapes" of the distortion (run "o2") , the energy radiated 
in £ = 5 can be as large as 21%. Similarly, runs corresponding to distorted Kerr black 
holes ("rO" to "r5") contain different mixtures of the different multipoles depending on 
the shape of the initial distortion and on the rotation rate of the black hole. As a trend, 
large rotation seems to yield a large i = 3 component: the maximum corresponds to run 
"r4", for which the i = 2 component carries 92. 1% of the energy and i = 3 carries ~ 7% 
of the energy. 

Three-dimensional, non-axisymmetric evolutions are urgently needed. A handful of 
non-axisymmetric simulations can be found in Ref. [24]. They find that the total energy 
radiated in a given multipole, Ef^, ranges from £'20 ~ 3 x lO^^M to £42 ~ 3 x lO^^M. 
Only modes with £ = 2, m = 0, 2 and £ = 4, m = 0, 2 should occur at linear order in the 
amplitude, all other modes being quadratic in the amplitude (for one of these "nonlinear" 
modes they find a completely negligible energy, £52 ~ lO^^^M). 

In fact, one of the most interesting features introduced by nonlinearities is mode- 
mode coupling. Ref. lEsIl studies this coupling in detail for nonrotating (Schwarzschild) 
black holes. Large- amplitude incident waves effectively add mass to the black hole, 
producing a small but visible "redshift" in the QNM frequencies. The waveform is also 
slightly shifted by nonlinearities. This nonlinear dephasing might be an important issue 
when applying matched-filtering techniques. On the other hand, when the amplitude 
A of incident waves exceeds a certain threshold, nonlinearities amplify the outgoing 
wave, which of course is good news for detection. Beyond the "nonlinearity threshold" 
the outgoing amplitude of the £-th multipole scales as A^/^. The reason is that, for 
ordinary spherical harmonics, ijirnt ~ ^(/t£)ra+tenTis with smaller I. Given a linear 
incident wave with £ = 2, an £ = 2fc mode arises from order k (and higher) nonlinear 
terms; this mode's amplitude scales as A^-f (^(A^+^). So, the quadratic coupling of an 
{£ = 2, m = 0) mode with itself produces an (£ = 4, m = 0) mode of order + G{A?) 
(and higher), in addition to an (£ = 2, m = 0) mode of order A -f ^(A^) (notice that 
£ = is ruled out for gravitational perturbations, being non-radiative). In addition, 
nonlinearities produce beating phenomena. In general, quadratic terms of an harmonic 
with frequency ft) contain the frequencies (0, 2o;), cubic terms contain (ft), 3oj), and 
so on. The selection rules for mode coupling in Kerr backgrounds should be more 



interesting: even in the linear regime, the Lense-Thirring effect couples axial and polar 
perturbations of a slowly rotating star according to Laporte's rule |26]. Mode couplings 
in the Kerr background can be very important for gravitational wave data analysis, and 
deserve further study. 

In conclusion, simulations suggest that (at least in situations with a certain degree of 
symmetry) high overtones are heavily suppressed. Investigating an extensive catalogue 
of initial data, Ref. |23] shows that we just don't know how energy is going to be par- 
titioned between different modes in a realistic merger event, unless we can predict with 
accuracy the shape of the initial data. Nonetheless, for data-analysis purposes it should 
be reasonable to assume that modes with £ = 2, m = 0, ±2 have the largest amplitude 
(the components with |m| =2 probably being dominant). The contribution from £ = 4, 
m = 0, ±2 could be smaller by about two orders of magnitude in energy. Modes with 
t — 3 may also be relevant, especially at large rotation rates. Nonlinearities could pro- 
duce a cascade of energy into higher-order modes, but we need more simulations with 
realistic initial data to assess the relevance of this effect for QNM detection. 

Head-on black hole collisions - 

In four dimensions, the total energy radiated b)^a freely-falling particle released from 
rest at infinity in the Schwarzschild spacetime [8, "221] is in surprising agreement with 
numerical simulations of head-on black hole collisions. An extensive comparison of 
numerical results with perturbation theory can be found in Refs. [28, 29, 30. 31. 32] . 

Fig. 14 and Sec. IV of Ref. ll30l] show that the perturbation-theory prediction for 
the radiated energy, is in excellent agreement with numerical 

results when we replace the particle's mass by the system's reduced mass. In fact, 
perturbation theory slightly overestimates the energy output predicted by numerical 
simulations, which (according to state-of-the-art simulations) is £ ~ 0.0013M ll33ll . 
In the so-called "particle-membrane" picture, this disagreement is compensated for 
multiplying the perturbative prediction by three fudge factors fsO*]: i) a factor iv^, 
accounting for the fact that in numerical simulations the infall starts at finite distance 
ro; ii) a factor F/, coming in because the black hole (unlike the falling particle) has 
finite size, and tidal deformations heat up the horizon; iii) a factor F^bs to account for 
reabsorption of the emitted gravitational waves by each black hole. If the collision starts 
at infinite separation, F^^ = 1; F/, has a weak dependence on ro, and tends to F/, ^ 0.86 
as ro — » oo; Fabs — 0.99 is even less relevant. 

If the total energy output from head-on black hole collisions is well understood, unfor- 
tunately the relative excitation of higher multipoles with respect to £ = 2 is not. Ref. |l32ll 
compares numerical calculations with the close-limit approximation (valid when the 
black holes are close enough to be considered as perturbations of a single Schwarzschild 
black hole) and with the particle-membrane approach. Fig. 12 of Ref. l32ll compares the 
radiated energy according to the close-limit approximation with results from the numeri- 
cal simulations. Extraction of the £ = 4 component from the simulations is very sensitive 
to numerical noise, especially when the black holes are very close and the £ = 4 compo- 
nent is more heavily suppressed. 

Numerical simulations have been extended to include unequal-mass black holes ll 34ll 
and boosted black holes [35, 36] . More recently, mesh refinement techniques were 
used to obtain more reliable estimates of the energy radiated in £ = 2 [331 . Even for 



highly boosted black holes, the agreement between perturbation theory (in the close- 
limit approximation) and numerical simulations is very good. Numerical studies of 
boosted black holes show that, fixing the initial separation between the holes, the energy 
radiated saturates at £ ~ O.OIM for very large values of the initial momentum of the 
holes (see eg. Fig. 2 of [36]). Ref. [35] combines Newtonian dynamics and numerical 
simulations of boosted black holes to suggest that the maximum energy emission could 
actually be much lower than this, E < 0.0016M. In general, these studies suggest that 
the role of the weak-field phase of the evolution is only to determine the momentum of 
the black holes when they start to interact nonlinearly, emitting most of the radiation in 
the ringdown phase. This expectation is confirmed by the comparison of perturbation 
theory and Post-Newtonian calculations the bremsstrahlung radiation emitted by 
a particle at separations larger than 4M (in Schwarzschild coordinates) contributes only 
~ 3% of the total energy, most of the radiation coming from the final (ringdown) phase. 



CONCLUSIONS 

Most numerical simulations so far consider nearly-axisymmetric situations, so they don't 
have much to say about the distribution of energy between modes with different values 
of (/,m) in a realistic black hole merger. This is a crucial issue for data analysis from 
ringdown waveforms [j|5|]. Recently there has been some remarkable progress in non- 
axisymmetric simulations of black hole mergers (see eg. [38] and references therein), 
but a discussion of these developments is beyond the scope of this paper. 

Some estimates of the energy radiated in the plunge phase by solar-mass binaries of 
rotating black holes are provided by the effective-one-body approach [39]. The error 
introduced by extrapolating to supermassive black holes should not be too large, since 
the energy balance of black hole binaries involves a single fundamental scale (their 
total mass M). It is useful to compare the effective-one-body estimate with earlier 
estimates of the energy radiated in the plunge plus ringdown phases, as provided by the 
Lazarus approach [W]. For aligned (anti-aligned) spins and specific angular momenta 
ji = j2 = 0.17 (ji = j2 = 0.25) the effective-one-body approach predicts an energy 
release ~ (0.6- 0.9) %M [~ (1 -3)%M]. The Lazarus approach finds ~ (1.7- 1.9) %M 
[~ (1.9 — 2.1)%M], respectively. In principle, the difference of these two estimates 
should bracket the most reasonable range of ringdown efficiencies. However, the overall 
error on both estimates is still too large to draw any definite conclusion. Numerical 
relativity should be able to provide a reliable answer to this problem within a very short 
time. 
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